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Abstract 

Heavy-baryon chiral perturbation theory (HBChPT) is applied to the asymmetry A^ in n+p — ► 
d + 7 at threshold, which arises due to the weak parity non-conserving interactions. Instead of 
appealing to Siegert's theorem, transition operators up to next-to-leading chiral order are derived 
and the corresponding amplitudes are evaluated with the Argonne vis wavefunctions. In addition 
to the impulse contribution, both parity-conserving and parity-non-conserving two-body one-pion- 
exchange diagrams appear up to this order. Our prediction for the asymmetry is A 7 = —0.10 h^ N , 
which is close to the Siegert's theorem based result, A^ ~ —0.11 h^ N . This illustrates that 
HBChPT is effectively applied to the parity-non-conserving physics. 
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There are still large discrepancies in the value of the hllr^, where is the leading order 

weak irNN coupling constant. For example, detection of the circular polarization of 1081 keV 
gamma rays from 18 F transition predicts \h^ N \ < 1.43 x 10~ 7 Measurement of anapole 
moment of 133 Cs predicts, on the other hand, much larger value, = (9.5 ± 0.22 ± 0.34) x 10~ 7 

|@]. The process 

n + p — > d + 7 (1) 

at threshold is free from the uncertainties of many-body (A > 3) systems, and thus a suitable 
source for the study of the parity non-conserving (PNC) effects in nuclear reactions. As a relevant 
observable for this purpose the asymmetry Ay is investigated, which is defined by the dependence 
of the cross section on the angle defined by the directions of the photon emission and the neutron 
polarization, W(9) oc 1 + Ay cos 6. The first calculation of the Ay was performed by Danilov || 
obtaining Ay = —0.08 hlljN- More elaborate attempts were made with realistic wavefunctions 



ill, 



Ay = -0.11 h^. (2) 



In more detail, Hamada-Johnston, Reid-soft-core and Tourreil-Sprung potentials are adopted to 
yield Ay = —0.109 h^ NN , —0.114 h^ NN and —0.107 h^ NN , respectively [[§. The available data are 
from the ILL experiment H, Ay LL = —(1.5 ± 4.8) x 10~ 8 . This data with the eq.(j2j) then imposes 
^nNN = (l-^ i 4-4) x 10~ 7 . At LANSCE [||, an experiment that aims at having 10~ 9 accuracy in 
Ay is under progress, which will sharpen the determination of greatly. 

The good convergence among the above various calculations is mostly due to Siegert's theorem 
|], 10 1 . It relates the major part of the El amplitude at low-energy to the one-body charge 
density whose amplitude can be estimated reliably without detailed informations of the nucleon- 
nucleon reactions. This indicates that a substantial departure from eq.(|2|) is unlikely. In the 
meantime, effective field theories (EFTs) have recently gained great successes in low-energy two- 
nucleon systems, which include the Solar proton fusion [11], the total radiative np capture cross 
section [12] and its spin observables [13, 14 1, the deuteron properties and low-energy nucleon-nucleon 
phase shifts ||, [Q| (FJ, |l|. Kaplan, Savage, Springer and Wise jOJ (KSSW) performed an EFT 
calculation of the process (||), using the so-called power-divergence subtraction scheme. They found 
rather surprising result, Ay = 0.17 h^ N , where the difference in overall sign is simply due to a 
mismatch in conventions (see Ref. [20, ||lj). Desplanques |2C] has analyzed KSSW's result in great 
detail. He showed that KSSW's result is - apart from the overall sign - exactly equivalent to the 
conventional result but with the zero-range approximation (ZRA) for the wavefunctions. The ZRA 
is responsible for the KSSW's overestimation of the asymmetry. 

In this paper, we will show that the asymmetry can be understood accurately by HBChPT, an 
EFT that has been thoroughly tested in low-energy nuclear physics. For this purpose we will go to 
next-to- leading-order (NLO) in Weinberg's power counting J2^], with the wavefunctions obtained 



by the Argonne vis potential [23]. So far this hybrid method has been found to be quite powerful. 
For example, the total cross section of the neutron thermal capture was found to agree to the 
experimental data perfectly with the theoretical error bar about 1% 12]. 



The leading-order PNC Lagrangian takes the form p, 24, |l 

72 



Cpnc = _MV e 3a6 iV t T V 6 iV, (3) 
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Figure 1: Generic diagrams that contribute to the asymmetry. The solid, dashed and curly lines 
are nucleons, pions and photons; and the "x" marks denote the insertion of the PNC vertex. 
Empty blobs stand for the sum of arbitrary number of iterations of the parity-conserving strong 
interactions, while the shaded blobs with (without) the "x" mark do the PNC (PC) two-nucleon 
irreducible currents. None of the crossing diagrams have been drawn. The shaded blob in (a) will 
be denoted as J, while that in (b) will be denoted by 5J. 



where e = +1 and h^j^ is the weak coupling constant that should be determined from the 
experimental value of the asymmetry. There are other PNC interactions at higher order, but 
our current poor understanding of the nuclear PNC effects does not make much sense to include 
them. Thus we will limit ourselves to the above leading order PNC term. Due to the smallness 
of h^lr N (~ 10~ 7 ), it is sufficient to consider only the contributions linear in I\1t N for the PNC 
amplitude. The factor h/^ NN appears either in the wavefunctions or in the currents, 

5{V f \J\*i) = (<5*/|Jl*i> + (*/|J]«5*t) + {*f\8j\*i), (4) 

where VE^ (^/) and J are the initial (final) wavefunctions and the electromagnetic currents, respec- 
tively. The "<5" marks stand for the first order perturbation with respect to the h^ N . Generic 
diagrams for the asymmetry are drawn in Fig. [l], where the first and the last terms in eq.(Q) corre- 
spond to Fig. (a) and Fig. (6), respectively. The second term, being symmetric to the first, is not 
shown in the figure. Note that the final state wavefunction in Fig. (a) has a parity-violating vertex. 

As mentioned, we follow Weinberg's power counting p^ ], where an irreducible diagram is 
counted as of order (Q/A x ) u , Q is the typical momentum scale and/or pion mass and A x ~ 4^/^ ~ 
rriM is the chiral scale. The chiral index v is given as 

v = 2L-2(C-l)-l + ^v h 

i 

where C and L are the numbers of the separate pieces and loops, respectively. A vertex indexed 
by the subscript i is characterized by v% = d{ + ^ + ej — 2, where dj, n.- L and e» are the number 
of derivatives/m^'s, nucleon lines and the external fields, respectively. For PC interactions, chiral 
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symmetry guarantees > [25]. Fig. (al) is LO with C = 2, while Fig. (a2, 6) are NLO with 
C = 1. However, the counting rules for the process at hand is tricky. For example, the PNC 
amplitudes due to Fig. (a) contain the two-nucleon reducible part, while the above counting rule is 
for irreducible diagrams. The unnatural smallness of the binding energy of the deuteron, B^, can 
also contaminate the counting rule, in the sense that a quantity proportional to positive powers of 
Bd becomes much smaller than what the counting rule implies. We find that the convergence of the 
chiral expansion is rather poor and we should go to at least up to NLO to have accurate results. 

To be specific in convention, we also write down the parity-conserving (PC) strong lagrangian 
explicitly, 

Cpc = -i^Na-Vit i T i N + ■■■, (5) 

with qa — +1.26, f n ~ +93 MeV and the ellipsis denotes terms not relevant in the discussion. 
Eqs.(||, H) then lead to the PNC one-pion-exchange (OPE) potential that is identical to the one 
given in @, |, |, 

M = ^HxW + 4)^(^), (0) 

where r = r\ — f 2 , r = \r\ and r = f/r. Note that when Siegert's theorem is used, the sign of the 
asymmetry is determined by the PNC potential. Thus the consistency of our convention to the 
conventional studies are guaranteed by the above PNC potential. ^ 
The one-body currents (Fig. (al)) read 



Jib = > — — i — -Oi x fe 7 H 

^— ; |_ 4 m 7v 4m n 



(7) 



where fis = fJ-p + Mn — 0.880, fly = Mp ~~ A*n — 4.706 and the ellipse denotes higher order terms. 
The two-body PC (Fig. (o2)) and PNC (Fig. (b)) currents read 

J2B = JWf ( l 2) WW ml + ? 

= -rSl^i x Ui • a 2 (m^ (r) - S®(f)) + S 12 (f)y 2 (r)} , (8) 

XT f ? h *NN 9A d {a x +a 2 )-q 

OJ 2 B — — I 7— \Tl-T% — TiTn 1= — - — ^rz, o ^ — 

J (2vr) 3 2V2 Udq ml + q 2 

h (l) 

= ~f ^^[n-T 2 -T^}(l + m 7r r)y (r)r-(a 1 + a 2 ) (9) 

did Q-m^r 

2/2(0 = r- — -z~yo(r), y (r) = — (10) 

or r or iirr 

Si 2 (f) = 3(?i • fa 2 ■ f - <j\ ■ a 2 . 



where 



and 



#5 KSSW [|19[ used PC Lagrangian with different sign, which causes different sign in the PNC potential and the 
asymmetry. 
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The PC currents, Jib and J2B, contribute to the asymmetry through the small parity-odd compo- 
nents of the wavefunctions induced by the PNC potential, as drawn Fig. (a). Here we recall that 
the PNC potential eq.(||) does not commute with either isospin or orbital angular momentum and 
induces parity-odd components (which are linearly proportional to h^ N ) in the spin-triplet wave- 
functions. On the other hand, the PNC currents, SJ2B1 contribute to the asymmetry by connecting 
the parity-even components. 

At threshold, there are only S-waves in the initial np states, ^{j 00 ( 1 S'o) and ^j* 1 ( 3 <Si), where 
^j^ J are the partial waves. We write the wavefunctions as 

= ^n s (r)CioXoo, 
v47rr 



S*° ; ;V1 = ^(<?i + ?2)-™,(rKioXu,, (11) 



where x(C) represents spinor(isospinor). We multiplied in front of the P— wavefunction to make 
vt(r) real at threshold. The deuteron wavefunction is the same with <J/ 011 but with the subscript "d" , 
instead of "t". The radial functions are normalized as linv^oo u s (r) = r — a s , Hindoo ut(r) = r — at 
and lim J ._> 00 Ud{r) = e _7d ' r , where a s ~ —23.7 fm (at — 5.42 fm) is the spin-singlet (spin-triplet) 
np scattering length and 7^ = \JBdrriN. The boundary conditions for the 3 Pi radial wavefunctions 
are lim r _ >0o v t (r) oc h^l N ^ and lim^oo v d (r) cx hy NN (l + ^ 



It is well-known (see, for example, |T^]) that the total cross section near threshold is predomi- 
nated by the isovector Ml transition, 1 Sq — > d, 

J\^°) = x\ Md [-K*p ~ ff n) x kM( l S )] xoo (12) 

with 

M( 1 5 ) = ^ [°°dru d (r) u s (r) + ■ ■ ■ = (1 + <5 2B ) (0.263 fm 2 ) (13) 

4771 jv JO 

where /c = wfc is the momentum carried out by the photon, the ellipsis denotes the two-body- 
current contributions and c>2B = (4.6 ±0.3) % denotes the ratio of the two-body currents compared 
to the one-body contribution. For the transition from the spin-triplet np state, there are both PC 
isoscalar Ml and PNC El contributions, 



(ttdlJM") = x\ Md -i(a P + *„)xfcM( 8 £i) 



3, 



XlJz 



6(y d \J)V jl l ) = x t 1Md [('? P + ^)E( 3 < S 1 )jxij z . (14) 

The isoscalar Ml transition is tiny (less than 0.1 %) in the total cross section and does not contribute 
to the asymmetry. Detailed analysis can be found in [13, 14 1 and we will neglect this isoscalar 
contribution hereafter. 
The asymmetry reads 

E( 3 ^) 

where the E( 3 Si) is the consequence of the PNC interaction and proportional to h^ N . The LO 
E( 3 Si) comes from the one-body (IB) contribution (Fig. (al)), while the two-body contributions 



A = -2^7T7TT> (15) 
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(Fig. (o2)) and (&)) are NLO, 



E( 3 ,Si 



il -kNN 



E — Eib + E2B 



(16) 



with 
Eib 
E2B 



2\/6mjv 
12x/6/2 



5.4 



v' d (r) \ u t {r) + 



Wt\Tj 

V2 



+ 



Vd{r) 



Jo 



mly {r) - 5 {3) (r) + 2y 2 (r) 
w d (r 



dr r 

dry (r)(l + m n r) [ u d (r) + 



u t {r) - V2w t (rf) + (d <-> t) 
w d (r) 



V2 



Ud(r) + 
u t {r) + 



V2 
w t (r) 
V2 



v t {r) + {d^t) 



(17) 



where v d j(r) is defined by v d j(r) = v d ^{r) /h^ NN and the u (d *-> i)" denotes the permutation 
between the subscript u d" and u t" . With Argonne wig wavefunctions we have E = (0.0428 — 
0.0302) fm 2 = 0.0136 fm 2 and consequently 



^ChPT 



-0.10 h 



(i) 

ttNN ■ 



(18) 



It might be worthwhile making a comparison of our result eq.(18) with the Siegert's theorem 
prediction, 



E 



Siegert 



10 

4V6 



dr r 



Vd{r) u t (r) + 



wt{r) 
V2 



(d t) 



0.0148 fm 2 



(19) 



where the numerical value is obtained with Argonne vis wavefunctions. Our result is close to the 
Siegert's theorem prediction, but the chiral convergence is rather slow, 



Eib + E2B 

-Siegert 
hi 



2.89- 1.97 = 0.92. 



(20) 



While a more study is needed to be definite, here let us present a plausible scenario for the bad 
convergence. One can see that the Siegert's theorem prediction, eq.(|T9|), is suppressed by the 
smallness of to ~ B d , which is smaller than its natural size, Q 2 /A x ~ ra\jm^ ~ 21 MeV, by about 
10 times. The suppression mechanism is, however, not manifest in our HBChPT results, eq.(0). 
To understand the consequences of this, let us expand Eib and E2B with respect to B d at threshold, 



E nB = A^ + B d A^ + 0(B 2 d ), (n=l, 2). 
Then what Siegert's theorem tells us is that 



i(0) 
L 1B 



+ A 



(o) 

2B 



0. 



(21) 



(22) 



In case both 4b and 4b 

are not zero, the above equation indicates that the ratio E2B/E1B 
becomes —1 at B d — * limit, which is not much far from our result, E2B/E1B — —2/3. In this 
respect the net result eq.(20) is the result of strong cancellation between the -B^-independent terms, 
A^Q- By the same token, the above scenario says that our PNC El amplitudes can be contaminated 
by a small deviation from eq.(^), when B d is quite small. However it should be understood clearly 
that the HBChPT results are to be improved systematically by taking the higher order. Thus, for 
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example, to include the two-body vector charge contributions which have been neglected in Siegert's 
theorem, we should rely on a systematic EFT like HBChPT. It is quite promising to observe that 
HBChPT up to NLO could already explain the A^/h^ N ratio within 10 % (compared to the 
Siegert's theorem prediction). We would like to make it clear that the situation is completely 
different in, for example, Ml transition amplitude, where the suppression factor u ~ can be 
factored out and we are left with the quantity which is non-zero even at = limit. In this case, 
a beautiful chiral convergence has been observed [Q |. 

So far, we have limited ourselves to the threshold limit, where only S np states are relevant. 
At thermal energy, there is tiny but non-zero contribution from the np ty 111 state. This 3 Pi state 
gives non-zero El transition amplitude that is independent of the h^ N and dependent on p, the 
direction of the relative np momentum. This small P— component gives rise to PC asymmetry and 
its magnitude was calculated in ]2(|. The PC asymmetry arises from the PC scalar h ■ (p x kj). 
While PNC asymmetry measures up-down asymmetry, PC asymmetry is the left-right asymmetry. 
The amount of the PC asymmetry is reported to be ~ 7 x 10~ 9 which can contaminate the exact 
measurement of Ay. This corresponds to the accuracy goal of the LANSCE of the order 10 -9 . 
Experimental considerations to discriminate these false signals were illustrated in ||. 

We thank B. Desplanques for his valuable discussions during this work. TSP is grateful to S.-i. 
Ando for his comments. Work of DPM and CHH is partially supported by the KOSEF grant No. 
1999-2-111-005-5, KRF grant 1999-015-DI0023, 2000-015-DP0072 and BK21 program. 
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